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FLUCTUATION THEORY OF CONNECTIVITIES FOR 
SUBCRITICAL RANDOM CLUSTER MODELS 
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Universitd di Bologna, Technion and Universite de Geneve 

We develop a fluctuation theory of connectivities for subcritical random 
cluster models. The theory is based on a comprehensive nonperturbative prob- 
abilistic description of long connected clusters in terms of essentially one- 
dimensional chains of irreducible objects. Statistics of local observables, for 
example, displacement, over such chains obey classical limit laws, and our 
construction leads to an effective random walk representation of percolation 
clusters. 

The results include a derivation of a sharp Ornstein-Zernike type asymp- 
totic formula for two point functions, a proof of analyticity and strict con- 
vexity of inverse correlation length and a proof of an invariance principle for 
connected clusters under diffusive scaling. 

In two dimensions duality considerations enable a reformulation of these re- 
sults for supercritical nearest-neighbor random cluster measures, in particular, 
for nearest-neighbor Potts models in the phase transition regime. Accordingly, 
we prove that in two dimensions Potts equilibrium crystal shapes are always 
analytic and strictly convex and that the interfaces between different phases 
are always diffusive. Thus, no roughening transition is possible in the whole 
regime where our results apply. 

Our results hold under an assumption of exponential decay of finite volume 
wired connectivities [assumption (1.2) below] in rectangular domains that is 
conjectured to hold in the whole subcritical regime; the latter is known to be 
true, in any dimensions, when q = 1, q = 2, and when q is sufficiently large. 
In two dimensions assumption (1.2) holds whenever there is an exponential 
decay of connectivities in the infinite volume measure. By duality, this includes 
all supercritical nearest-neighbor Potts models with positive surface tension 
between ordered phases. 
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1. Introduction and results. We recall the definition of random clus- 
ter measures P in Section 1.1 below. Let x £ be a distant point. The 
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event {0 x} means that x is connected to or, equivalently, that the 
common cluster Cq^ x of the origin and x is well defined. In the subcrit- 
ical situation the probability P(0 <-> x) decays exponentially with \x\. The 
direction-dependent decay exponent £ is called inverse correlation length. 
Our main structural result in this paper is a probabilistic description of 
connected clusters Co, 2 in terms of one-dimensional "sausages" of objects, 
which we call irreducible clusters; see (1.4) below. Under the conditional 
measure P(- | <-»• x), the sizes of these irreducible clusters have exponen- 
tially decaying tails. This is the original Ornstein-Zernike picture. It was 
rigorously justified for several particular models either in a perturbative 
regime in, for example, [5, 19, 20] or for on-axis directions in, for exam- 
ple, [7, 12]. In the latter works the arguments heavily relied on particular 
microscopic structure of irreducible clusters and on independence between 
different irreducible components. The main contribution of [8] is a deriva- 
tion of robust coarse-graining procedures which fits in with essentially any 
microscopic notion of irreducibility one prefers to choose. The issue of de- 
pendence between irreducible components has been resolved for all subcrit- 
ical finite-range Ising models in [9], where the corresponding statistics have 
been studied in terms of one-dimensional thermodynamics of Ruelle's type. 
This is a generic situation. Below we extend the theory of [8, 9] to a large 
class of finite-range subcritical random cluster measures. The theory ap- 
plies in any dimension. For two-dimensional nearest-neighbor models dual- 
ity arguments enable an extension for supercritical interfaces. In particular, 
our theory leads to a comprehensive fluctuation analysis of interfaces for 
any two-dimensional nearest-neighbor Potts model in the phase segregation 
regime. 

1.1. Random- cluster measures. Let J = (Jx) X £Z d De a collection of non- 
negative real numbers satisfying J x = J_ x ; we suppose that there exists 
R > such that J x = as soon as \x\ > R (\x\ denotes Euclidean norm), 
and that 3 J x > when |x| = 1. Let us denote by £j the set of all unordered 
pairs (x,y) of sites of Z d such that J y - X > 0; we call bonds the elements of 
£j, and endpoints the two elements of a bond. 

Let f2 = {0, be the set of all bond configurations, and let J- denote 
the corresponding product cr-algebra; we say that a bond e is open, respec- 
tively closed, in the configuration u € O if u(e) = 1, respectively u>(e) = 0. 
A configuration of u G O is identified with the graph with vertex-set Z d 
and edge-set {e € £j :oj(e) = 1}. The maximal connected components of this 



3 This second assumption is purely a matter of convenience; the proofs still apply, 
modulo straightforward modifications, if this condition is removed, provided that the graph 
(Z d ,£j) remains connected. 
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graph are called open clusters of the configuration (in particular, each iso- 
lated site is considered to be an open cluster). We write x <-> y for the event 
that x and y belong to the same open cluster, A <-> B if there exist x £ A 
and y S B such that x <-> y, and x <-> oo for the event that the open cluster 
containing x has infinite cardinality. 

If A is a finite subset of Z d , A <g Z d , we denote by £ j the set of all bonds 
with both endpoints in A, and by 7a the cr-algebra generated by the sets 
{w(e):ee£j\S^}. 

The random cluster measures on Z d , with coupling constants J and param- 
eters j3>0 and q> 1, are all the probability measures P on (0, T) satisfying 
the following DLR equations: \/A 6 J^" and A d Z rf , 

P(A | Ta)^) = f>l ;JAq (A), for P-a.e. r? G Q, 
where P^-j/jg is the probability measure on (O,^ 7 ) given by 

0, otherwise, 

with p e d = l-e- 2f3J e, Q v A = {u}en:u(e) = r)(e),Vee£j\££}, and *rt A (w) 
the number of open clusters of cj intersecting A. 

We write P^.j ^ when 77 = 1 (wired boundary conditions), and P^.j ^ 
when 77 = (free boundary conditions), and shall generally omit J,/3,q from 
the notation. 

The measures P A and P^ are defined for each finite AcZ* 1 and can be 

cA 

viewed as distributions on the set {0, 1} J . In a completely similar fashion 
one constructs measures P^ and P™ on {0, l} 7 for any finite subgraph (set 
of bonds) 7 of (Z d ,£j). 

A property of random cluster measures that is repeatedly used in the 
present work is the FKG property: Let / and g be two increasing functions 
w.r.t. the natural partial order on O; then ¥(fg) > P(/)P(g) for any random 
cluster measure P. 

In particular, it follows that, for any increasing sequence A / Z d , the 
two sequences of measures P^ and P A converge to random cluster measures 
P w and P f . Moreover, the following stochastic ordering of measures obtains, 
P =3! P =3! P w for any random cluster measure P. 

It can also be shown that PJ^ = Pj ^ (and therefore, there is a unique 
random cluster measure), for all but at most countably many values of (3 [15]. 

In particular, the following definition, 

/? c ld =l f inf{/3:P} i/3i(? (0^oo)>0}, 

does not depend on the particular choice * = w or * = f (since the event 
<-> 00 is increasing). The set {f3:(3<(3l} is called a subcritical regime. 
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It can be shown that there is a unique random cluster measure for each 
P < Pi [15]. 

From now on, we suppose that < and denote the corresponding 
(unique) random cluster measure by Fj,/3, g , or simply P when no confusion 
is possible. 

1.1.1. Inverse correlation length. The central quantity in our investiga- 
tion is the connectivity function P(0 <-> x). Associated to the latter is the 
inverse correlation length: for x € R d , 

f (ac) d = - lim ~ logP(0 <-> [kx]), 

where [y] £ Z rf is the component- wise integer part of y 6 R d . 

In the subcritical regime, the connectivity function decays with \x\. It is 
natural to ask whether this decay is actually exponential. We thus introduce 
another notion of critical point: 

/? c 2d =l f sup{/3: min £(n)>o). 

Obviously, (5^ < fi\ , but it is actually believed that these two critical points 
always coincide. 

Conjecture 1. For any d>l, q>l, and J as above, 

This result is known to hold when q = 1 [1], q = 2 [2] or g is large 
enough [17] (in any dimensions). 

Let us assume now that ft < 01 . The function £ is clearly positively homo- 
geneous; moreover, FKG inequalities imply that £ is a convex function on 
R d , and £ is obviously finite. It is therefore an equivalent norm on R d . Two 
convex bodies play an important role in the sequel: the equi-decay set (or, 
equivalently, Frank diagram) and the Wulff shape 4 Kg, which are defined 
by 5 

(1.1) V^ f {xER d :£(x)<l}, K^ f p| {teR d :(t,n) d <m}- 

Ties' 1 - 1 

Notice that these two sets are polar, 

Ug = \x£R d :max(t,x) d <l\, Kg = It £ R d : max(t,x) d < ll. 



4 This terminology will find its justification when stating our results in the 2D Potts 
model. 

5 (x,y)d denotes the scalar product of x, y € M. d . 
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Pairs (x,t) G R d x <9K^ are called dual if (t,x)d = £(x). The set (or, 
equivalently, U^) encodes geometrically all the information on the function 
£ (in particular, the latter can easily be reconstructed from the set). 

1.1.2. Duality in dimension 2. We restrict here our attention to the two- 
dimensional nearest-neighbor case: J x = lrui =1 i; in this case we simply write 
£ instead of £j . 

A specific feature will allow us to reinterpret some of our results as per- 
taining to the g-state Potts model in the phase coexistence regime: duality. 

Let Z 2 '* = {(x*, y*) : x* - \ , y* - \ G Z} be the dual lattice, and denote by 
£* the nearest-neighbor bonds of Z 2 '* (the dual bonds). Let A = {— n, . . . , re} 2 

and A* = {— n — i,.. . ,n + ^} 2 C Z 2 '*. Random cluster measures on £1* = f 
{0, 1} £ * are defined exactly as before. 

There is a natural mapping from $7^0*: given a configuration well, 
the configuration uj* G W is obtained by setting, for each dual bond e* G £* , 
uj(e*) = 1 — u>(e), where e G £ is the bond dual to e* (i.e., seen as unit-length 
line segments in M 2 , e is the bond intersecting e*). 

Let p = f 1 — e -2 ^, and define p* as the solution to the equation p*/(l — 
p*) = q(l - p)/p; let (5* be such that p* = 1 — e~ 2/3 * . Then the following 
duality relation holds [13]: 

PW^)= P A*;/3*>*)- 

Of course, letting n — > oo, one then also obtains the corresponding statement 
for infinite-volume measures: q {A) = P^* ^(^4*), for all cylindrical events 

A and A* = {u*:uj£A}. 

This duality relates super- and subcritical models. In particular, this will 
allow us to derive, from our results obtained in the subcritical regime, results 
about the supercritical regime. Thus, the inverse correlation length £ of the 
subcritical model coincides with the surface tension (see, e.g., [17, 18]) r of 
the dual supercritical model. Much more than that: since connected clusters 
of the subcritical model represent interfaces in the supercritical dual model, 
our results lead to a comprehensive fluctuation theory of microscopic inter- 
faces which applies to all nearest-neighbor two-dimensional random cluster 
measures (at q > 1), in particular, it applies to all nearest-neighbor two- 
dimensional Potts models in the phase segregation regime. 

1.1.3. The basic assumption. Although the natural region of validity of 
our results should be the whole region {/?:/?< /3 2 }, where £ > 0, it will be 
necessary to have an a priori stronger property. We therefore introduce yet 
another critical point f3 c , defined as the supremum over all values of (3 such 
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that the following holds: Let Ayv = {— N, . . . , N} ; then there exist vq and 
v\ > such that, for any N, 

(1.2) P^^^O-S^A^Ooe-^. 

Obviously, /3 C < (3%, and it is not difficult to show that $ c > 0, for all d > 1, 
q > 1, and J as above. Again, the following is expected: 

Conjecture 2. For any d>l, q>l, and 3 as above, 

h = vl- 

This result has been proved in all the cases mentioned above, where it 
is known that fi\ = (5^: q = 1, q = 2 [6], and q sufficiently large [17], in any 
dimension. Moreover, Conjecture 2 has been verified for any q > 1 in two 
dimensions [3]. In the Appendix we give a simple proof of the latter assertion 
which was kindly explained to us by Reda Messikh. 

Our basic assumption throughout this paper is that (3 < (3 C . 

1.2. Effective one- dimensional structure of long connected clusters. From 
the technical point of view, we encode long connected clusters Co,x as one- 
dimensional strings of letters from a countable alphabet of irreducible clus- 
ters. The main effort is to verify that the resulting one-dimensional structure 
falls in the framework of Ruelle's thermodynamic formalism for full shifts 
with Holder continuous potentials. Once such identification is accomplished, 
statistical properties of various local observables comply with classical limit 
laws. In fact, the only local observable we consider here is the displacement 
over an irreducible cluster. This leads to an effective random walk represen- 
tation of Co,x and provides adequate tools for proving all the main results 
of the paper. However, we would like to point out that our method provides 
much more flexibility, and that there are many other local and quasi-local 
observables of interest. 

Thus, before stating precise results, let us give a somewhat informal de- 
scription of the above mentioned probabilistic structure of clusters Co,x in 
terms of a string of irreducible clusters, which lies in the heart of our ap- 
proach: With each t G OK^, we associate three families 2T t b ,Xt and I t of 
respectively backward irreducible, irreducible and forward irreducible clus- 
ters. These families are precisely defined in Section 2.10; see also Figure 1. 
An important geometric feature is the cone confinement property, which we 
proceed to describe. 

1.2.1. Geometry of irreducible clusters. For each t G , define the cone 

Y = {y£Z d :(t,y) d >K 2 ay)}- 

Then t-irreducible clusters satisfy the following Y-cone confinement condi- 
tion (see Figure 1): For any 7 b G T t b , there exists b G 7 b , for any j f G 2^, 
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there exists f G 7 f and, finally, for any 7 G li, there exist b,f € 7 such that 

(1.3) 7 b Cb-Y,7 f Cf + y and 7 C (f + Y) n (b - Y). 

Notice that two irreducible clusters 7 and 7' with b(7) = f (7') could be 
patched together. We use 7 II 7' for the corresponding concatenation. 

In (1.4), (1.6) and (1.8) below we state some of the main properties which 
hold uniformly in x G Z d and t G 3K^. Namely, there exist two constants 
Kl = «l(J ; P, q) > and K2 = «2(J, /3, (?) G (0, 1) such that: 

1.2.2. Decomposition of Co jX . Up to probabilities of the order 

exp{-(t,x) d - m\x\} 
clusters Co, x can be represented concatenation, 

(1.4) c , a; = 7 b U7in---n7Arn7 f , 

with M > 1, 7 b G J f b , 7 f G ^ and 71, . . . , 7A r G T t . 

Of course, such a statement is meaningful only for x's satisfying £(x) < 
(t,x) d + /ci|x|. 

In Section 3.1 we construct percolation events T b , Ti, . . . ,T_\f and T f , so 
that the probability of each fixed realization of (1.4) is precisely 

(1.5) P(r b r x • • • ivr f ) = P ( r b )p(r0^^^p. 

1.2.3. Weights of boundary conditions. Boundary clusters 7 b and 7 f 
cannot be too large: Uniformly in y G Y , 

* * 

(1.6) J2 F(T b )<e- {t ' y)d - KM and ]T P(r f ) < e -<*.v)rf-«a|l#l , 

where the summation is restricted to backward and forward irreducible clus- 
ters satisfying b(7 b ) = and, accordingly, f (j f ) = 0. 
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1.2.4. Statistics of the effective random walk. For the purpose of this 
paper, it is convenient to think about the decomposition (1.4) in terms of 
an effective random walk with boundary conditions 7 b and 7 f and steps 

V( 71 ),...,V( 1M ), 

where V(j) = b — f is the displacement along the irreducible cluster 7 G It] 
see Figure 2. In this notation, 

V{j) = V{ 11 ) + --- + V{ 1M ) 

is the total displacement along the AT-step effective random walk. 

For each couple 7 b and 7 f of boundary conditions, there exist two func- 
tions g{- I 7 b ,7 f ) and tpt(- | 7 b ,7 f ) acting on strings of irreducible clusters 
from It, such that the fraction on the right-hand side of (1.5) can be rep- 
resented in the following way: Let 71 , . . . , 7^ be irreducible clusters and let 
Ti, . . . , Tj^f be the corresponding percolation events: 

P^Ti-iyi*) (tyh))d 
p(r b )p(r f ) 

(1.7) = exp{V>*(7i, . • • ,7V I 7 b ,7 f ) H ^MlAf I 7 b ,7 f )} 

x 5(7i 5 •••,7A/' I 7 b 5 7 f )- 
The potentials ipt satisfy 

e V , t(7iv,7AT|7 b ,7 f )+Kl|V(7i)| < ^ 

7i62t 

uniformly in t G 9K^, in 7 b G 2^, 7 f G 2"^, G N and in 72, ... , 7at G X t . 
In the Bernoulli percolation case (q = 1) [8] both ^ f and 5 depend only on 
the first element of the string. This corresponds to effective random walks 
with independent increments which, by (1.8), satisfy Cramer's condition. 
Accordingly, the fluctuation analysis of connected clusters boils down to an 
application of classical local limit results of probability theory. For general 
random cluster models, both the potentials ipt and the functions g depend on 
all the variables involved. Nevertheless, in Sections 3.3 and 3.4 we show that 
they possess appropriate regularity properties, which enable an equally clas- 
sical local limit analysis along the lines of one-dimensional thermodynamics 
of full shifts over countable alphabets. 

1.3. Main results. The geometric description (1.4) of Co, 2 in terms of 
one-dimensional chains of irreducible clusters and the regularity properties 
of potentials in the induced effective random walk representation (1.7) (see 
Sections 2 and 3 for precise statements) provide a unified framework for all 
the results below. 
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Fig. 2. Effective random walk representation. 

1.3.1. Ornstein-Zernike behavior of connectivities. The principal quan- 
tity of interest in this work is the connectivity function, 

P(0^x), 

and, more specifically, its asymptotic behavior as \x\ — > oo. Our first result, 
which extends earlier results in the case of self-avoiding walks [16] , q = 1 [8] 
and q = 2 [9], is that these asymptotics are of the Ornstein-Zernike type. 

Theorem A. Let (3 < /3 C . Then 

P (° " x ) = | g ^-ij/2 exp(-g(n,)M)(l + o(l)), 

uniformly as \x\ — > oo. T/ie functions ^ and £ are positive, locally analytic 
functions on S d_1 , and = £c/|x|. 

As a corollary to the previous theorem, we obtain the following sharp 
asymptotics for probabilities of exits from sets. For simplicity, we shall for- 
mulate such result for cubes A at and under the assumption that J is invariant 
under all reflection symmetries of Z rf , however, a generalization to generic 

piece- wise smooth domains Vjv = NV and general J is straightforward. 6 

Corollary 1.1. Assume that J is invariant under all reflections ofL d . 
Let (3 < [3 C . Then there exists a constant ip, such that 

(1.9) P(0 Z d \ A N ) = 2d^ N ^ ei) {l + o(l)) . 

1.3.2. Geometry of Wulff shapes and equi-decay profiles. As mentioned 
above, there is a deep relationship between the geometry of the sets and 

and the inverse correlation length £. As a by-product of our analysis, we 
obtain the following information on the latter. 



Here generic means that the minimum of £ is attained on a discrete set of points of 
dV, and dV is smooth around these points. 
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Theorem B. Let (3 < (3 C . Then Kg has a locally analytic, strictly convex 
boundary. Moreover, the Gaussian curvature xp of Kg is uniformly positive, 

d-i 

where {X/3i(t)} are the principal curvatures o/<9Kg at t. By duality, dUg is 
also locally analytic and strictly convex. 

1.3.3. Invariance principle for long connected clusters. The local limit 
result stated in Theorem A can be extended to a full invariance principle 
for the common cluster Co :Xn of and x n G Z rf as x n — > oo. The magnitude 
of fluctuation depends on the asymptotic direction of x n . Let x G S d_1 and 
let t G (9Kg be the dual point. Consider a sequence of vertices x n = [nx\ 

and the corresponding sequence of conditional measures ¥ x n 
x n ). By our basic cluster decomposition result, Co,x n has the form (1.4) 
up to P X;ri -probabilities of order e~ Kin . Let 7 = {0, Uq, . . . , uj\f, x n } be the 
trajectory of the corresponding effective random walk, and let C n [^\ be 
the linear interpolation through the vertices of 7. Let 7i x be the (d — 1)- 
dimensional hyperplane orthogonal to x. Alternatively, 7i x is the tangent 
space of <9Kg at t. By the cone confinement property (1.3), the intersection 
number 

#(£ n [ 7 ]n(rLnxJ +H X )) = 1 

for every r G [0, 1]. Accordingly, there is a natural parametrization of C n [y] 
as a function $ n : [0, 1] 1 — > Ti x . Namely, 

*„(r) =proj Wx (£ ri [7] n (r[nx\ +H X )). 

We shall define the diffusive scaling of Co jX „ in terms of the diffusive scaling 
of $ n . For r G [0, 1], set 

n (r) = -=S„(r). 

Such definition of diffusive scaling of Co l3 ; n is justified since by (1.6) and 
(1.8) the Hausdorff distance between Co,x n an( l £-n[l] satisfies 

(1.10) lim P x , n (d H (C ,^,/:„)>^logn) = 0, 

n — ^00 

for some K = K{(5) large enough. 

Theorem C. Let (3 < (3 C . Then {4> n (')} weakly converges under {¥ x , n } 
to the distribution of 

(y/xJJB^-), y/X04-iBd-x(-)), 
where B\, B^-i are independent standard Brownian bridges on [0,1]. 
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1.3.4. Two dimensional Potts models. In the special case of two-dimen- 
sional models with nearest-neighbor interactions J x = l{i x i = i}, duality al- 
lows us to reinterpret the preceding results as results about interfaces in the 
supercritical random cluster model, as explained in Section 1.1.2. In partic- 
ular, when q is an integer, this provides a detailed description of interfaces 
in the g-state Potts model in the phase coexistence regime. Let us briefly 
recall the definition of this model. Let q > 2 be an integer, and A <s Z 2 . 
The nearest- neighbor g-states Potts measure in A with boundary condition 
a € {1, ... , g} z , at inverse temperature j3, is the probability measure on the 
product u-algebra of {1, ... , q} z2 given by 

f (Z%)- 1 \{ if a = a off A, 

r ' C ki. a ) = { {i,j}nA/0,|i-j'|=l 

1. 0, otherwise. 

Since the inverse correlation length of the subcritical model coincides with 
the surface tension of the supercritical model, the set really corresponds 
to the Wulff shape associated to the dual model (see [4] for a review on Wulff 
construction). Therefore, observing that, as already mentioned, (3 C is known 
to coincide with /3 2 , we see that we can restate Theorem B as follows. 

Theorem D. Potts models in two dimensions do not have roughening 
transition. In other words, for any q £ N and at any temperature (3 at which 
the surface tension £ of the two-dimensional nearest-neighbor q-state Potts 
model is positive, the corresponding equilibrium crystal (Wulff) shape 
is strictly convex and has an analytic boundary with strictly positive 

curvature xp- 

This result was previously known only in the cases q = 1 and q = 2. 

Let us now consider the 2D g-state Potts model in Ajv = {— N, . . . ,-/V} 2 , 
with n-Dobrushin boundary condition, n G S 1 : cfj = 1 if (ft, i)d > 0, and 
Oi = 2 if (n,i)d < 0. For definiteness, let us assume that (n,e 2 ) d > l/y/2, 
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Fig. 3. The interface of the 2D nearest-neighbor three state Potts model. 
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and write for the corresponding measure. Under such a boundary con- 
dition, we can define unambiguously the interface, for any configuration a 
having positive probability under H^: Among all connected components of 
all dual edges separating disagreeing spins, only one has infinite cardinal- 
ity; the intersection of the latter with the box A defines the interface (see 
Figure 3). It turns out (see the proof of the following theorem) that the 
interface as defined here is a subset of an open cluster in the subcritical dual 
of the random cluster model to which it is associated through the standard 
coupling [21], and thus can be approximated by the same effective random 
walk. Denoting as before by $tv its linear interpolation, and by 4>n the dif- 
fusive scaling of the latter, we arrive at the following theorem (notice that 
this is not an immediate corollary of Theorem C because the corresponding 
random cluster model is defined in a vertical strip, and one must therefore 
take care of the boundary effects). 

Theorem E. Consider the q-states nearest-neighbor Potts model with 
n-Dobrushin boundary condition in the box An, at inverse temperature (3 
such that the surface tension is positive. Then {4>n(')} weakly converges 
under to the distribution of 

tyxW)}> 

where B is the standard Brownian bridge on [0, 1] and X/3 is the curvature 
of dK.£ computed at the point t £ dK.^ dual to n. 

The case q = 2 has already been treated in [14] . 

Furthermore, the theory we develop also leads to a nonperturbative mi- 
croscopic justification of the Wulff construction for two-dimensional Potts 
models on the DKS-level, and also implies that the inverse correlation length 
is an analytic function of (3. However, we relegate the corresponding discus- 
sion to a future paper. 

1.4. Structure of the paper. As it is indicated in Section 1.2 our main 
new result in this work is a robust stochastic-geometric description of long 
connected clusters Co,x, which is valid for all subcritical random cluster 
measures satisfying assumption (1.2). 

In Section 2, a finite scale renormalization analysis of the typical geometry 
of connected clusters implies that Co,x has, up to negligible probabilities, 
an essentially one-dimensional structure, which can be conveniently visual- 
ized as a chain of irreducible clusters. The crucial model-oriented input for 
an adaptation of the general coarse-graining techniques developed in [8, 9] 
is a mixing estimate on exponential decay of connectivities under various 
boundary conditions. This estimate is derived in Section 2.1 as a direct con- 
sequence of our assumption (1.2). The rest of Section 2 is devoted then to 
the construction of the irreducible decomposition (1.4) of the cluster Co a:- 
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The second step of the proof is a refined local limit analysis of chains 
of irreducible clusters. Section 3 is devoted to a representation of P(0 <-> x) 
in terms of thermodynamics of Ruelle operators with uniformly Lipschitz 
continuous potentials for full shifts over countable alphabets (of irreducible 
clusters). The main representation formula (3.3) is derived in Section 3.1. A 
mapping to the thermodynamic formalism is explained in Section 3.3. The 
key result is Theorem 3.1 of Section 3.4 which asserts that assumption (1.2) 
implies appropriate regularity properties of potentials. 

Once Theorem 3.1 is established, the general theory of [9] applies and all 
the results announced in Section 1.3 follow in a relatively straightforward 
fashion. In Section 4 we briefly explain the corresponding proofs along the 
lines of [9, 10, 14]. 

Finally, in the Appendix we give a simple proof that in two dimensions 
assumption (1.2) holds whenever the inverse correlation length is positive. 
(This proof was explained to us by Reda Messikh and we publish it here 
with his kind permission.) 

Remark on constants. Positive constants c\ , C2, ■ ■ ■ are used in the inter- 
mediate computations and are updated with each section. Constants k%,K2 
and u±,U2,. ■■ appear in our principal estimates and are fixed throughout the 
paper. Finally, finite scale renormalization procedures we employ depend on 
three parameters K, r and 8 as it is described in Section 2 below. 

2. Geometry of typical clusters. 

2.1. Notation and basic decay estimate. The geometry of the problem is 
conveniently spelled out in terms of the equi-decay set defined in (1.1). 
Let us fix a number r > and a finite scale K > 0, both of which will be later 
chosen large enough, depending on (3 and on the cone-opening parameter 5 
which we shall define in Section 2.6 below. For any y G set 

B K (y) d ^ (y + K-XJs)nZ d and B K (y) d ^ B K+rlogK (y). 
For y = 0, we shall use the shorthand notation 

B/<= f B^(0) and, accordingly, = f B^-(O). 

Furthermore, for any subset iCZ 1 *, set 

A K d = U B r i ogK (y). 

y&A 

We shall also use the notation 

d R A d ^{y£Z d \A:d(y,A)<R} 

for the .R-outer boundary of A, where R denotes the range of the process 
(i.e., the length of the longest possible bond). Given a subset A C % d and 
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two points x G A and y G A U <9#^4, we shall use {x y} to denote the 
event that x and y are connected by an open path X:x *— > y, such that all 
the vertices of A, with the possible exception of the terminal point y itself, 
belong to A. 

There are two parameters K and r (which enter the picture through the 
definitions of B^, B/^ and Ak above) to play with. Let us say that a certain 
quantity is of order o^-(l) if for every power m > one can find r > and 
a scale Kq , such that for every K > Kq this quantity is bounded above by 
1/K m . 

The next proposition is the only place where we rely on our basic esti- 
mate (1.2) (beyond just assuming exponential decay of connectivities). 

Proposition 2.1. For any subset A C B^ and every vertex y G 8rBk, 

(2.1) P w f(0 ^y)< e -^(1 + Oi ,(l)). 

Proof. The proof is rather standard (see Figure 4): By the FKG prop- 
erty of P, we can decompose P^(0 y) as follows: 

P I A . (° " v) = ¥ A K (o ^ v\ 9rA <- d R A K ) 

+ FJ k (0 ^y;d R A^d R A K ) 
< P ^( y)FJ KXA (d R A ~ d R A K ) 
+ P^(0^y). 

By assumption (1.2), 

FJ KXA (d R A~d R A K ) = o K (l). 

On the other hand, 

2.2. Skeletons. In this subsection we shall develop a coarse-grained pic- 
ture of the common cluster Co tX of and igZf Fix a scale -fT and a number 
r as in the beginning of Section 2.1. 

Given a realization C of Co a, we are going to construct a tree 7~(C), with 
vertices xo = 0, xi, . . . , a;jv(T) usm § the following procedure. 

Step 1. Set x = 0, C = B K (x ), i = 1. 

Step 2. If there is at least one vertex y G d R C such that 

y^d R B K (y)\C, 



16 M. CAMPANINO, D. IOFFE AND Y. VELENIK 




Fig. 4. The event {0 <- — > y}. The dashed curve is the boundary OrAk- 




Fig. 5. The tree skeleton associated to Co, x . 



then choose y* to be the minimal (in the lexicographic order) such vertex 
and go to Step 3. Otherwise stop the procedure. 

Step 3. Set Xi = y*. Update C — > C U B K (xi) and i — > i + 1. Go to Step 

2. 

This procedure produces the vertices of a tree T(C). Clearly, 
[_jB/^(xi) is connected, and C C [J &2K(%i)- 

i i 

An example of a tree skeleton of a cluster is depicted on Figure 5. 
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We now give a rule to construct the edges of the tree: By definition, the 
vertex xi is connected to the vertex of 

{xj,j <i:xi£ d R B K (xj)}, 

which has the smallest index. We define the probability of such a tree T as 

P(T) djf £ r( Co , x = c), 

where the sum is over all clusters containing both and x compatible with 
the above procedure (in the sense that the tree T can be obtained from the 
cluster using the procedure). 

By Proposition 2.1, we have the following bound on the probability of 
observing a particular tree T: For each i = 0, 1, . . . , N(T), set 

^ = B^(xi)\ljBx(^). 

j<i 

By our construction of T, 

P(T) < F( Xi ^ d R B K (xi) for i = 0, 1, . . . , N(T)) 

N(T) 

(2.2) < J] PJ^( Xl ^d R B K ( Xl )) 

<eM-K(l-o K (l))(N(T) + l)), 

where we have used (2.1) on the last step. 

We are now going to relabel the vertices of a tree T in the following way: 
T can be conveniently split into a trunk t = (to, ■ • ■ ,tjv(t))) an d a family of 
disjoint branches 55 = (b , b 2 , . . .); T = (t, 55). The trunk is defined as the 
path in T connecting to = xo to t;v = xp, where the latter vertex is defined 
to be the vertex X{ of T with the smallest index i, having the property that 
x £ B2_R-(xj). Each branch b l is then a connected sub-tree rooted at some 
vertex of the trunk. The symbols N(t) and 2V(Q5) = J2i N (b l ) 

are reserved 

for the total number of vertices of t and 55 respectively. Of course, N{T) = 
N(t) + iV(55) in this notation. 

The probability of a trunk t is defined as 

P(t) d ^ f ^P(C 0i:c = C), 

where the compatibility C ~ t means that t is a trunk compatible with this 
construction, and similarly for the leaves. The main idea of the renormal- 
ization approach is that on large enough finite scales K typical trees have 
a simple geometry: Most of the increments along the trunk t point into the 
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direction of the target vertex x, whereas the branches 23 are very short and 
sparse. Thus, our estimates on probabilities of trees follow the scheme 

(2.3) P((t,»))=P(t)P(® |t). 

Our bounds on P(QS | t) are crude: One pays at least e~ K Q~°KW) for each 
additional vertex of T and, as we shall see in Sections 2.3 and 2.4, at large 
enough scales K such a price beats the entropy of different tree patterns. 
The main point is to control the forward geometry of typical trunks t, which 
we shall do with the help of the surcharge function introduced in Section 2.5. 

In all the estimates below probabilities of various tree patterns will be 
tested against an a priori bound 

(2.4) P(0^x)x e - ?(x) . 

2.3. Typical trunks do not contain too many vertices. The first observa- 
tion is that the trunk t of Cq )X cannot contain too many vertices. 

Lemma 2.1. There exist r, a scale Kq, and two constants c\ and C2 > 
such that 

V(t)>ciM|0^^ <e~ c ^, 

uniformly in x and K > Kq . 

Proof. The number of trunks of m vertices is bounded above by 
(c 3 (d, J R)J^~ 1 ) m , and therefore, using (2.2), 

Hm>c4)< £ „ P (- ra *(i- ,(i)-M£i^>) 

m>c\\x\/ K 

and the conclusion follows by (2.4). □ 

2.4. Typical trees do not contain too many branches. 

Lemma 2.2. Let k > 0. There exist r, Kq(k), and C4 = C/i(k) > such 
that 

f(n(<B) >Kjjr\°^ x ) <e" C4K|:r| , 
uniformly in x and K > Kq . 

Proof. By Lemma 2.1, we can assume that N(t) < ci^. We shall call 
such trunks admissible. Given an admissible trunk, the number of trees with 
N(*B) = b> is bounded above by (see Lemma 2.3 in [8]) 

exp (c$b (log (-*) + log if 
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Therefore, for any admissible trunk t with N(t) < ci -feL, one has 

p(tf(!B)>«M|l)< £ exp(-«(l-c 5 !2i(l^t^ 

fo^** Ki | x 1 1 a. 

< e -«N/2 ; 
provided -ftT is sufficiently large. □ 

2.5. Surcharge function. On large enough renormalization scales -ftT, most 
of the increments along typical trunks t point in the direction of x. More 
precisely, let t G dK.^ be a dual vector such that (t,x)d = £(x). Then the 
probability of a trunk t= (to, . . • ,tj\r) is quantitatively measured using the 
surcharge functional 

N 
Z=l 

where the surcharge function St : M d i— > K + is defined via 

s*(y) d =£(y) - (t,y) d . 

By construction, Sj(x) = 0. 
Now, for all y G <9j?B^, 

if > £(y) -c 6 r log K = {t,y) d + s t (y) - c 6 r log K. 

Since by construction t; G dp^xih-i) for every / = 1, ...,JV(t), and since 
x G B2A-(tjv), we infer, setting n = iV(t), that 



K > ]T^£(tj - ^-i) ~ cernlogK = (t, t N ) d +s t (t) - c 6 rnlogK 
l 

> £(x) +St(t) - c 7 rnlogK, 

where on the last step we have tacitly assumed that ti» K. Reasoning as 
in (2.2), we arrive to the following upper bound: 

(2.5) P(t) < exp (-fa) - a t (t) + c 8 ^|^|x|) , 

for every admissible trunk t. 

We are ready to state our crucial surcharge inequality, which enables a 
control of the forward geometry of typical trunks of Co tX : 

Lemma 2.3. Let e > 0. There exist r and Kq{e) such that 
¥{5 t {i)>2e\x\^^x)<e-^ x) - £ \ x \ 

uniformly in x G vectors t G <9Kg dual to x, and renormalization scales 
K>K . 
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Proof. By Lemma 2.1, we can restrict attention to admissible trunks. 
The number of the latter is bounded above by 

\x\ logK^ 



exp[c 9 (d, R)- 

Therefore, in view of (2.5), 

P(*t(t) > 2e|x|,JV(t) < ci\x\/K,0*->x) 

<exp(-\x\(2e- (c 8 r + c 9 )^^^ 

and the conclusion follows once K is large enough. □ 

2.6. Cone points of trunks. Let 5 E (0, g). For any i £ <9K^, we define the 
forward cone by 

Y>(t) d ^{xeZ d : St (x)<5Z(x)}, 

and the backward cone by Y^(t) = f — Y^(t). 

Given a trunk t= (toj • • • jtjv(t))) we say that t& is a (t, <5)-forward cone 
point of t if 

{t k+1 ,...,t N(l) } Ct k +Y 5 >(t), 
while tfc is a (i, <5)-backward cone point of t if 

{to, . . .,t fe _!} c t fc + Y<(t) = t fe - y 5 >(t). 

tfc is then said to be a (t,5)-cone point of t if it is both a (i, <5)-forward cone 
point and a (t, <5)-backward cone point of t. 

If t contains points that are not forward cone points, we define 

11 = min{j : tj is not a (t, <5)-forward cone point of t}, 
r> = min{j > l x : tj - t h $ Y>(t)}, 

12 = min{j > r\ : tj is not a (t, <5)-forward cone point of t}, 
r> =mm{j>l 2 :tj-ti 2 ^Y>(t)}, 



and, similarly, if t contains points that are not backward cone points, 
If = max{j :tj is not a (t, <5)-backward cone point of t}, 
rf = max{j < h : tj - t h (£ Y s K (t)}, 

I2 = max{j < r\ \Xj is not a (t, 5)-backward cone point of t}, 
= max{ j < l 2 : tj - t h £ Y 5 < (t)}, 
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We then say that tj is a (t, <5)-marked point of t if 

jey{l>,...,r>-l}u\J{r< + l,...,l<}, 

k k 

where V denotes the disjoint union. 

The next lemma, the proof of which is the same as the proof of Lemma 2.2 
in [9], controls the number of marked points, #™ < f rk (i), in terms of the sur- 
charge function. 

Lemma 2.4. Let 5 G (0, |) be fixed. The surcharge cost St(t) is controlled 
in terms of the number of marked points as 

uniformly in x, dual directions t £ OK^ and large enough values of r and K . 

In view of Lemma 2.3, we infer that, there exists c\\ such that, for every 
e>0, 

(2.6) P(#27 k (t)>e^|0~x) <exp(- Cll eS\x\), 

as soon as K > K (e,5), for all dual pairs x,t. In particular, since for any 
trunk t we have a trivial a priori bound 

(2-7) N(t) > Cl2 M 

inequality (2.6) implies that most of the vertices of typical trunks on large 
enough scales are cone points. We can formulate the latter assertion as 
follows: There exists a constant C13 > such that for every e > one can 
find r and a finite scale Kq = Kq(£,5) such that 

(2.8) p(#™ ark (t) > eiV(t) \0^x)< e - c i3 £ M, 
for all dual pairs x,t and scales K > Kq. 

2.7. Cone points of trees. We now want to extend the previous estimate 
on trunks to an estimate taking also into account the branches. In order to 
do so, we have to slightly enlarge the opening of our forward and backward 
cones. Given a tree T = tU 23, let us say that a vertex tj of its trunk t = 
(to, . . . , ij, ■ ■ ■ , tj\r(t)) is a cone point of T if 

(2.9) T C (t,- + Y<(t)) U (t, + Y>(t)) = (t 3 - Y>(t)) U (t, + Y>(t)). 

Evidently, each cone point of T is also a cone point of t. The converse is, in 
general, not true: For tj E t, let us denote by b(tj) the branch of 25 which is 
rooted at tj. By Lemma 2.2, b(tj) = for the majority of i's. Some branches 



22 M. CAMPANINO, D. IOFFE AND Y. VELENIK 

are, nevertheless, nonempty, so let us say that a cone point tj of t is blocked 
if (2.9) does not hold. Equivalently, a cone point tj of t is blocked if there 
exists tj E t such that 

(2.10) b(tf) 2 (ij + Y< s (t)) U (t 3 +Y>(t)). 

We claim that the number #^ ocked (t, «8) of all blocked cone points of t is, 
with overwhelming probability, small relative to the total number of all cone 
points of t. Precisely, the following uniform (in x, dual directions t £ 9K^, 
large enough scales K and admissible trunks t) bound holds: 

Lemma 2.5. Let an admissible trunk t satisfy #%f Tk (t) <eN(t). Then 
there exists C14 > such that 

(2.11) P(#^ ockcd (t,<B) > eiV(t) 1 1) < e - Cl4E|x| . 
Proof. As in (2.2), 

P(b(to) = 60, • • • , b(tjv) = b N I t) < exp^-(iT - o K (l))J^N(bi)J , 

for each particular realization {bj} of {b(t,)}. For every i = 0, . . . , iV(t), define 
Xj to be the number of cone points of t which are blocked by b(tj), 

X % = #{j : (2.10) holds}. 

Clearly, #^ ocked (t, «B) < Ei^i- On the other hand, in view of the (cone) 
confinement property of cone points of T, there exists C13 > such that 

N(b t ) > c 15 Xi. 

Therefore, exactly as in the proof of Lemma 2.2, we conclude 
¥(X 1 >x u ...,X N >x N \t)<ex^-^-K^2x l \ 

and (2.11) follows. □ 

2.8. Typical tree skeletons of Co tX . We summarize the results of 
Sections 2.3-2.7 as follows: Let x, t and 5 be fixed. Given a tree skele- 
ton T = (t, 03) of Co tX , let us denote by t* = (t^, . . . , U N * ) the collection of 
all (t, <5)-cone points of T. Here N* = N*(T) stands for the total number of 
these points. 

Theorem 2.1. For every 5e (0,5), there exists r, a scale Kq and two 
positive numbers v% and u±, such that 



<-> x \ < e 



-U4\X\ 



(2.12) F (^ N *(T)<^ 
uniformly in x, dual directions t £ <9K^ and finite scales K > Kq. 
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2.9. Cone points of Co,,. We say that a vertex y G Co,, is a (t,<5)-cone 
point of the latter if 

C ,, C {y + Y<(t)} U{y + Y>(t)} = {y- Y>(t)} U {y + Y>(t)}. 

Notice the triple cone opening in the definition above. 

If we choose 5 to be too small, then it might happen that Y^(t) does 
not contain axis directions, and, consequently, the above definition might be 
meaningless. Since, however, £ is an equivalent norm, it is not hard to see 
that one can choose 5 G (0, |) in such a fashion that for every t G dK.^ there 
exists an axis direction e*j such that either efj or — ei belongs to the interior 
of Ygg(t). Let us fix such a value of 5 for the rest of the paper. 

Any cone point t of T is a reasonable candidate for a cone point of Co,, 
itself. In fact, since by construction the Hausdorff distance between T and 
Co, s does not exceed some c±qK, and, furthermore, away from the origin 
the cone Y^(t) lies in the interior of the double cone Y^(t), there exists a 
constant cn, such that tj is a cone point of Co,, whenever 

Bcrffe) n Co,, C {t,- - Y>(t)} U {tj + Y>(t)}. 

By the finite energy property of P, the probability of the latter event is 
bounded below by a positive number p = p(K, 5) regardless of the percolation 
configuration on the complement Z rf \ B Cl7 #-(tj). Therefore, for a fixed tree 
T which contains a subset t* = (t*, . . . , t* M ) of cone points which in addition 
satisfy 

B ci 7 tf(C+i) n B cirK (t* m ) = 0, 

the P(- | 7")-conditional distribution of the number of cone points of Co,, 
stochastically dominates binomial distribution Bin(M,p). 

As before, let t* = , . . . , ti N „ ) be the collection of all cone points of 
T . By Theorem 2.1, we can restrict attention to the case when N*(T) > 
v^\x\jK. By construction, all the increments tj, — tj, G ^w^) whenever Z > k. 
Consequently, 

l-i 



e(ti, - ti j > (*, ti, - 1» j d = £ (*, ti m+1 - tj m ) d 
> (i - 1) J2 e(t tm+1 - o > a - fc)(i - 



m=k 



Accordingly, 



B Cl7 i^ (Ui ) n B ci 7 i^ (tifc ) = whenever l-k> 2 ( Cl7 + _ 

l — o 

Let #j5 nc (Co,,) be the number of all (t,6)-cone points of Co,,. We have 
proved the following: 
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Theorem 2.2. There exists 5 G (0, |) and two positive numbers v§ and 
v§ such that 

(2.13) H#tT( C o,x) <vb\x\\0++x)< e -*W 
uniformly in x and in the corresponding dual directions t G . 

Theorem 2.2 paves the way for an irreducible decomposition of Co,z, which 
we proceed to develop. Apart from (2.13), none of the estimates of Sections 
2.3-2.9 is ever used in the sequel, and we shall feel free to reset the values 
of numerical constants c\-cyj. 

2.10. Irreducible clusters. Let 5 be fixed as in Theorem 2.2, and let 
t G 5K^. A finite -R-connected subgraph 7 C Z d is called backward (t,5)- 
ir reducible, 7 G Z t b , or just backward irreducible, if there is no ambiguity 
about t, if there exists a vertex b = b(7) G 7, such that 

(2.14) 1 Cb + Y<(t) = b-Y>(t), 

and, in addition, b is the only cone point of 7. 

Similarly, 7 is called forward irreducible, 7 GT{, if there exists a vertex 
f = f (7) G 7 satisfying 

(2.15) lQf + Y>(t), 

and, in addition, f is the only cone point of 7. 

Finally, 7 is called irreducible, 7 G It, if there exist f 7^ b G 7 such that 
both (2.14) and (2.15) are satisfied and, in addition, f and b are the only 
cone points of 7. 

In view of Theorem 2.2, we are entitled to ignore clusters Co,x which 
contain less than two (t, <5)-cone points (assuming that \x\ is sufficiently 
large, of course). The remaining clusters Co tX admit the unambiguous irre- 
ducible decomposition (1.4), where 7 b is backward irreducible, 7 f is forward 
irreducible, whereas the remaining clusters 71,..., 7/v are irreducible. The 
number M of irreducible components depends on Co,x- The notation 7117' 
means that the two clusters are compatible in the sense that b(7) = f(7'). 

3. Probabilistic description of typical clusters. In this section we show 
that typical long percolation clusters have a one-dimensional structure which 
complies with thermodynamical formalism of full shifts over countable al- 
phabets. 
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3.1. The representation formula. Let t and 5 be fixed so that (2.13) 
holds. With each (t, <5)-irreducible cluster 7 G I t , we associate the percolation 
event r = Wr\N = WN ', which is defined as follows: Let £w(l) be the edge 
set of 7. Then, 

W = {All edges of £w{l) are open}. 
Accordingly, let us define the set of edges 

(3.1) £ N (j) = {(u, v) : u G 7} \ OMt) U f "(f) U £ + (b)), 
where 

r(f) d ^{(u ) f):uef-y 3 >( i )} 

and 

£+(b) = {(b,v):veb + Y>(t)}, 

and, of course, f = f (7) and b = b(7). Then, 

N = {All edges of £n{i) are closed}. 

Notice that we refrain from closing edges from £~(f) and £ + (b) in (3.1) in 
order to be compatible with the decomposition (1.4) of Co,x- 

The definition of percolation events r b = W b AA b and F = W f AA f , which 
correspond to backward and, respectively, forward irreducible clusters, is 
completely similar. 

In this notation (1.4) gives rise [up to factors of order 1 + o(e _I ' 6 ' a; ')] to 
the following representation formula: 

(3.2) no ~x)=J2 E F(r b r x ---r n r f ), 

7 b 907 f 3x n 7l>— i7n 

where the last sum is over all strings of irreducible compatible clusters 7 = 
(71, ... , 7 n ) (and with respect to associated percolation events Ti, . . . , F n ). 

It is convenient to consider all the clusters in question modulo their Z rf - 
shifts. In this way any collection 7 b ,7 = (71, . . . ,7n)>7 f of, respectively, back- 
ward, string of irreducible, and forward clusters can be patched together. 
Introducing an additional notation 

y(7)= f b(7)-f( 7 ), 

for the displacement along an irreducible cluster 7, and the half-spaces 

ny^{yeZ d :(t,y) d >0} 

and 

H^- = {yeZ d :(t,x) d >(t,y) d }, 
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we rewrite (3.2) as 

(3.3) = J2 E E E p(r b )P(r f )e ( *' :c -^ +J/)d 



xE E 



where the shifts of backward and forward clusters are normalized in such a 
way that b(7 b ) =f( 7 f ) = 0. 

By Theorem 2.2, the above equality holds up to factors of order 1 + 
o(e-^l x l). 

3.2. Decomposition of probabilities. We shall decompose the probabili- 
ties in (3.3) as follows: First of all set 

r'Ti • • • r n r f = w b AA b WiM • • • w n M n w f M f d = w b N b wN, 

where W = Wi • • • W f and M = Mi ■ ■ ■M i . 

At this point we shall rely on a specific decoupling property of FK mea- 
sures 7 : namely, 

n 

P(W b W I M b M) = <bn 71 n...ii 7n ii,y (W h W) = P^ b (W b )pf f (W f ) E[P^(W fc ). 

i 

In the same fashion, 

p(r b )p(r f ) = p^ b (w b )Py(w f )iP(AA b )p(AA f ). 

Therefore, 

p(r b r x • • • r n r f ) " . . . w{M h Mi ■ ■ ■ M n M f ) 



p(r b )p(r f ) 11 i' 1 kJ 



P(AT b )P(AT f ) ' 

Finally, we decompose the above fraction as 

F(M b Mi • • -M n M f ) 
F(M b )¥(M { ) 

m = F(Afb| ^)^ f) w i * • ■ ■ W> • ■ - nv. i a"). 



7 Notice that the rest of our arguments depends only on positive association, finite 
energy and assumption (1.2). 
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Alternatively, we could have decomposed them as 

F(Af h M 1 ---M n Af f ) 
P(M b )F(M f ) 

(3.5) 

P(AA f | N n ■ ■ ■ M\Af b ) mt Kr ikf KrKrh , nfkrik/tt> 



P(AA f ) 



\N n | M n -i ■ ■ -NiN ) ■ ■ ■ P(M | AH. 



A comparison between (3.4) and (3.5) enables to formulate results in a par- 
ticularly convenient symmetric form. 

3.3. Reduction to one- dimensional thermodynamics. With clusters 7 b 
and 7 f fixed, define the function 

(3.6) g{iu ■ ■ ■ ,in) = g(iu ■ ■ ■ ,in 1 7 ,7 ) = pTj^b) ' 

and a potential 

(3.7) e <M^<-^) = e ^(7i,-,7n|7 f ) 4jf e^vM)^^^^ | ^ • • -Af n Af f ). 

Both the function g(-) and the potential tpt(') ac t on any string 7 of length n 
of (It) irreducible clusters for n = 1, 2, . . . . For any n 6 N, define a measure 
I 7 b >7 f ) on ^e set H n of strings of ra-irreducible clusters 

/4(71, • • • ,7n I 7 b ,7 f ) = e ^(Ti^--T»)+^^--T»)+-+^^) 9 ( 7l , . . . , 7n ). 

In this notation the last sum in (3.3) is just 

(3.8) fA(V(li) + • • • + V(ln) = z - y I 7 b , 7 f )- 

It is easy to see that Theorem 2.2 implies that there exists uj > 0, such that 

(3.9) ]T P(r b )e ( *' u)[i < e~ U7M and ]T P(r f )e (t ' u)d < e~ v ^ u \ 

uniformly in u^7L d . Consequently, the main contribution to the right-hand 
side of (3.3) should come from the terms (3.8), with z — y being close to 
x and with the number of steps n being close to optimal. The local limit 
analysis of measures /^(- | 7 b ,7 f ) boils down to a study of the analytic 
spectral properties of the Ruelle type operator 

(3.10) £*/(7)= E e^)/( 7 ,7> 

7 irreducible 

In this respect we shall be able to rely on a general theory developed in 
[9, 10] as soon as we check appropriate regularity properties of g(-) and 
tpt(')- This is the content of the next section. 
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3.4. Regularity properties of g(-) andtpt(')- For two strings of irreducible 
clusters 7 = (71, . . . , 7 n ) and A = (Ai , . . . , A m ) , define 

i(7, A) = f min{i : A,, 7^ 7^} A n A m. 

Theorem 3.1. There exist 6 < 1 and two positive constants v% and vg, 
such that the following inequalities hold uniformly in t G <9K.£, t-irreducible 
clusters 7, strings of t-irreducible clusters 7, A, t-backward irreducible clus- 
ters 7 b and pairs of t- forward irreducible clusters 7 f ; 7 . 

(3.11) ^8<5(7l7 b ,7 f )<-- 
Furthermore, 

(3.12) |5(7,7 I 7 b ,7 f ) "5(7,A| 7 b ,7 f )l < ^' A) , 
as well as 

(3.13) |^(7,7 I 7 f ) " Ml, A I 7 f )l < ^' A) - 

Proof. Given a set of bonds A, define 

AA(A)= f|{«(e)=0}. 

Thanks to the cone confinement property of irreducible clusters, all three 
claims (3.11), (3.12) and (3.13) of Theorem 3.1 are straightforward conse- 
quences of the following 

Proposition 3.1. Assume that (3 < 0%. Then, there exist a constant c\, 
such that 

uniformly in (finite) sets of bond A, B and C . 



vid(e,BAC*) 



• 



It is enough to consider C = B V D. Then, for every bond e, the FKG 
property of P implies 

< P(w(e) = I M{B V D)) - P(u(e) = | jV(JB)). 

However, a repeated application of the FKG property of P implies 

P(w(e) = I M(B)) > P(w(e) = 0; e ^ D \ M(B)) 

> P(w(e) =0] ey^ D;M(B))F(e j4 D) 

> P(cu(e) = j AA(5 V £>))(! - i/ e _1/ld(e ' D) ), 
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where we have also used exponential decay of connectivities on the last step. 
As a result, 

< P(w(e) = | M(B V D)) - P(w(e) = | M(B)) 
< z. e- lyid(e ' D) P(w(e) = | M{B V D)). 

Thus, 

n^e)=V\N{B)) _ Vld{e , D) 

-p{uj(e)=0\Af(BVD))- 

and (3.14) follows for any singleton A = {e}. 

The general case reduces to the singleton one. Let us number the bonds 
of A as A = {e±, . . . , e r }. After that set A k = {e±, . . . , e^}. Clearly, 

F(Af(A)\Af(B)) _ r TT 1 P(o;(efc + i) = | N{A k V B)) 
F(M(A) | M{B V D)) AA P( w ( efc+1 ) = | M(A k VBVfl))' 

and (3.14) follows. □ 

4. Implications. We are ready now to explain all our main results. 

4.1. Ornstein-Zernike formula for two point functions. Theorem 3.1 im- 
plies that the local limit analysis of (3.8) falls into the framework of the 
general local limit theory developed in Sections 4 and 5 of [9]. In particular, 
only terms with 

(4.1) \(z-y)-x\<g:\x\ 

really contribute to the right-hand side of (3.3). Furthermore, there exist a 
function F on the set of all backward and forward irreducible clusters, which 
is bounded above and below, 

-<F(-)<c 2 , 

C2 

for some C2 > 1, and a locally analytic positive function ^ defined on a 
neighborhood of ri x in § d ~\ so that 

E f&(V(n) + ■■■ + V{ln) = z - y | 7 b , 7 f ) 

(4.2) " . 

= paSp e-^+^-^F( 7 b )F( 7 f )(l + o(l)). 

Since —£(z — y) + (t, z — y)d < 0, (3.9) implies that the range of y and z 
which essentially contributes to the right-hand side of (3.3) can be further 
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shrunk to \y\, \ z — x\ < C3 log \x\ for some C3 large enough. For such y and z, 
however, 



i{z -y)-(t,z- y) d = i{z -y)- £0) - (*, (z-y)- x) 
uniformly in \x\ — > 00. As a result, the assertion of Theorem A follows with 



E nr f ) e (M)d F(7 f )) 



where we have relied on lattice symmetries of the model, and in particular, 
on the central symmetry employed in the definition of backward and forward 
irreducible clusters, on the last step. 

4.2. Exit from boxes: Proof of Corollary 1.1. The proof of Corollary 1.1 
is very similar to the derivation of Ornstein-Zernike asymptotics, so we only 
sketch the argument. 

First of all, observe that Kg inherits all the symmetries of J. Therefore, 
its strict convexity implies that the minimum of £(x) on {x€M d :|x|oo = l} 
is attained exactly when x = ±ej, i = 1, . . . , d. 

Next, one has the following obvious lower bound on the probability of 
exiting the box A^r: let x N = (N + l)ei, then 

P(0 Z d \ A N ) > P(0 x N ) > CN-^-^ 2 e- N ^ Sl \ 

for some C = C(@,J,d) > 0. Let t G <9Kg be dual to x N (notice that t is a 
multiple of e*i). Using the previous lower bound, one obtains that 

P(0 Z d \ A N ) = 2d¥(0 Y> 5 {t) \ Ajv)(1 + o(l)), 

so the problem is reduced to the analysis of the probability in the right-hand 
side. 

As before, up to factors of order 1 + o{e~ mN ), clusters connecting to 
the exterior of the box Atv, in the general direction specified by the forward 
cone Y^(t), admit a decomposition of the type (3.3), 

= E E E r(r b )F{T f )e (t ' xN - z+y)d 

y,z&A N 7 b 9 _ 3/7 f.7^( 7 )>(^ lj . r iv_ z ) d 

x V V P(r b Ti---r n r f ) (tz _ v)d 

p(r b )p(r f ) 
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where the horizontal range of a cluster C is defined as 
TZ(C) = sup(e*i,n) d - inf (ex,u) d . 

It is convenient to write z = (x^ — z±, z±), with z\ E Z and z± E Rea- 
soning as in Section 4.1, we see that one can restrict attention to points y 
and z such that \y\ < C3logiV and z\ < C3logiV, for some sufficiently large 
constant C3. In this case, writing z = (N + l,z±), we obtain that 

Y.V t n{V{ll) + --- + V{ln)=Z-y\ 1 \ 1 i ) 

n 

= ^^^^^ F ^ b )F^){l + 

Rough estimates show now that the contribution of terms with \z±\ > N l l 2+Ti , 
r] > 0, is negligible. For the remaining terms, we have 



E 

\z ± \<N 1 / 2 +V 



gKgg) r -e(z)+j(e 1 )N 



= V e~ N ^ d2 ^i^ 2 (l + o(l)) 

N (d-l)/2 e l 1 "^ 1 ^ 

= *(ei) J d d ^u e- (lt ' d2 « l =i lt)d/2 (l + o(l)), 
which completes the proof. 

4.3. Geometry of Wulff shapes and equi-decay profiles. By the relation 
between £ and K^, 

£0) = max(i, 

the Wulff shape Kg can be alternatively defined as the closure of the con- 
vergence set 

|t:^P(0^x)e ( *' x)d <oo|. 

Consequently (see Section 3.3 in [9]), given any t E dK.^, there exists e > 
such that the portion of the boundary inside the e-ball t + B £ around 
t can be described as the level set 

<9K € n (t + B e ) = {t + s : p(£ t ,s) = 1}, 



•32 



M. CAMPANINO, D. IOFFE AND Y. VELENIK 



where p(£t, s ) is the spectral radius of the operator 

Of course, £t jS is a perturbation of the operator £t defined in (3.10). In view 
of Theorem 3.1, the operator £ t satisfies all the assumptions of general The- 
orem 4.1 in [9]. In particular, p(£t) is a nondegenerate eigenvalue of £f and 
the Fredholm spectrum of £t lies strictly inside the C-ball of radius /o(£j). 
Local analyticity of 9Kg around t follows then from analytic perturbation 
theory of nondegenerate eigenvalues. Positivity of Gaussian curvature follows 
from the conditional variance argument. We refer to Sections 5.4 and 5.5 of 
[9] for the corresponding details. 

4.4. Invariance principle for long connected clusters. In view of the de- 
composition (1.4) and the regularity properties of the potentials as stated 
in Theorem 3.1, the proof of Section 3 of [14] just goes through. 

4.5. Nearest neighbor Potts models in the phase transition regime. Using 
the standard coupling of the Potts model with n Dobrushin boundary condi- 
tion and the associated random-cluster measure [21], it is easy to check that 
in any pair of configurations {uj, a) contributing to the joint measure one has 
the following property: All bonds dual to those belonging to the interface 
of a are closed in u (indeed, they join sites with spins taking different val- 
ues) . By duality, this means that all the bonds belonging to the interface are 
open in u*. Therefore, the distribution of the bonds in the interface is actu- 
ally stochastically dominated by the distribution of the bonds in the cluster 
C x n x n connecting the dual sites x^ d = (— N — |, — i) and x^, = f (N + ^, |) 
under the measure 

Pi(- 1 

at inverse temperature (3* (i.e., in the subcritical regime). The theorem would 
therefore follow if we can extend Theorem C to this finite-volume situation. 
As in [14], the only required change to the proof is to redefine the two 
extremal irreducible pieces 7 b and 7 f : one introduces two new extremal 

pieces 7 b = f 7 b II 71 II • • • II 7L(io g 7V)2J and 7 f d = TV-UiogTV) 2 ] U • • • U 7AT U 7 f • 
The size of 7 b and 7 f can be assumed to be at most O(logiV) 4 and therefore 
does not affect the limiting behavior, but is still large enough to keep the 
boundary effects under control. We refer to the detailed discussion in [14] 
(dealing with the case of the Ising model) . 



CONNECTIVITIES OF SUBCRITICAL RANDOM CLUSTER MODELS 33 

APPENDIX: BASIC ASSUMPTION IN TWO DIMENSIONS 

In two dimensions our assumption (1.2) holds once the inverse correlation 
length is positive. 

Proposition A.l. Let {J,(3,q} be fixed and assume that there is a 
unique {J , (3 , q} -infinite volume measure IP which, in addition, has exponen- 
tial decay of connectivities, 

P(O^x) <e" Cl|:r| . 
Then assumption (1.2) is satisfied. 

The above proposition has been established in a somewhat more general 
context in [3] . Here we sketch a simplified proof which was kindly explained 
to us by Reda Messikh (see also [11]) with a reference to the original idea 
by David Barbato. In order to facilitate the notation, we shall consider 
the nearest-neighbor case, so that the assumption on positivity of inverse 
correlation length is equivalent to the assumption on positivity of surface 
tension in the dual model. An adjustment to the general finite-range case is 
straightforward. 

Step 1. Define 

V(x) =I{ x ^dA N }- 

Then, J2 x eA N 7 ?( :E ) * s J us t the size of the boundary cluster. We claim that, 
for every 5 > 0, 



-oo. 



(A.l) limsuplpW-L J2 V(x)>s) 

The super-exponential bound (A.l) enables to restrict attention to arbitrary 
small boundary cluster densities. In fact, a stronger claim is true: 

(A.2) i imsup _Lpw f J_ £ V ( X ) >5}<0, 

for every 5 > fixed. (A.2) follows from the Hoeffding inequality: First of 
all, 

hm PX M (0~M M ) = 0, 

by the unicity of the infinite volume measure P. Thus, for every 5 > 0, one 
can choose M = M(5) < oo, such that 

(A.3) E Am(^ E Wm})<^ 
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Without loss of generality, we may assume that M(5) divides N, N = (2K- 
1)M. Thus, 

A N = V (2Mx + A M ). 

x£A K 



Since 

1 ^ / x ^ 1 ^ ( 1 x , 

L {y~d(2Mx+k M )} I , 



J_ V r,(x) < J- V ( J- V 



zSAjy 1 1 x€A K V JU 1 y€2Mx+A M / 

it remains to use (A. 3) and positive association of random cluster measures. 

Step 2. Next one constructs a minimal section type splitting of the 
event {0 «-> <9Atv/ 2 }- In the sequel we shall use the notation 

A l N d = Atv \ Ajv-j-i- 
Fix a small 5 > 0. By the preceding step, we may assume that 

5N 2 



(A.4) J2 V{x)< 



xeA N 



N/2 



In particular, at least one square layer (below | • | is the maximum norm) 

L l N = f {x : \x\ = N — I}, 1 = 1,..., N/2, has less that 5N sites connected to 
the outer boundary OAn. It is convenient to decouple the corresponding 
events from the event {0 <-> dAj^/ 2 } as follows: Consider random variables 

r a 1 " 1 

u(x) = \l, ifx^8A N , 
{ 0, otherwise, 



and set 



n(i)= Yl w(x). 

|x|=JV-J 



Of course, u>(x) < rj(x), and 



= min{Z : < 6N} e (l, - 



N 



for any bond configuration satisfying (A.4). Define the event A[ = {I* = I}. 
Then, ignoring configurations which violate (A.4), we arrive to the following 
decomposition: 

N/2-1 

{0^dA N/2 }= \/ {0~dA N/2 ;Ai}- 
1=1 
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Finally, note that Ai depends only on the percolation configuration on bonds 
(x, y) with both end-points belonging to A l N . For each 1 = 1,..., N/2 — 1, the 
latter set of bonds is disjoint from the set of bonds with both end-points in 

Aat/2- 

Step 3. Assume A\ and fix a bond configuration <pi on A l N and a bond 
configuration ipi on Ajv-z-i- Since there are less than 5N sites of the layer 
Ln_i which are connected (in <j>{) to dAjy, the finite energy condition implies 
that 

n N (<Pi^i) < e c " SN FJ N (cp l ^i;dA N/2 t4 dA N ) 
for some universal constant c 2 . It follows that, for each I, 

FJ N (0 <- dA N/2 ;A l ) < e c * SN F% N (0 - dA N/2 ; A N/2 *4 dA N ;A t ). 
However, 

FJJO <-> dA N/2 ;A N/2 8A N ) < P(0 8A N/2 ) 
by a standard FKG decoupling argument. 
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